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Abstract 

We study the cluster mass function and its evolution in different models with Dark 
Energy arising from a self-interacting scalar field, with Ratra-Peebles and SUGRA 
potentials. Computations are based on a Press & Schechter approximation. The 
mass functions we obtain are compared with results holding for open models or 
models with Dark Energy due to a cosmological constant. Evolution results, in 
some Dark Energy models, closely approach open models. 
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1 Introduction 



One of the main puzzles in modern cosmology is the nature of Dark Energy 
(DE), whose presence seems to be required by SNIa data (see, e.g., Perlmutter 
et al 1999, Riess et al 1998). A joint analysis of CMB and LSS observations 
(see, e.g., Percival et al. 2002, Efstathiou et al 2002) also favor a fiat Universe 
with a matter density parameter fi^ ~ 0.3, mostly due to CDM and with a 
minor contribution of baryons {Qbh"^ ~ 0.02; h is the Hubble constant in units 
of 100 km/s/Mpc; in this paper we shall take h = 0.7 and, unless differently 
specified, Qm = 0.3, anywhere). The residual energy content of the world, in 
the present epoch, should not be observable in the number-of-particle repre- 
sentation. 
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One of the most appealing possibilities is that such dark component arises 
from a self-interacting scalar field. With in the wide set of interaction potential 
suggested, a particular relevance is kept by Ratra-Peebles (1988, RP hereafter; 
see also Wetterich 1995) and SUGRA (Brax & Martin 1999, 2000) expressions: 

V{4>) = K^+^/r , V{(t>) = (A^+70°) exp(«:0V2) . (1) 

Here A is an energy scale, currently set in the range 10^-10^'^ GeV, relevant for 
fundamental interaction physics; potentials depend also on the exponent a; 
fixing A and a, the DE density parameter Vt£,E is automatically fixed; in this 
work we preferred to use as free parameters A and VLde', ^^ SUGRA potentials, 
n = SttG {G: gravitational constant). 

In this work we try to determine some effects on galaxy clusters and their 
evolution, caused by replacing a simple cosmological constant with DE due to 
a scalar field self-interacting according to the potentials (1). 

The technique used to study non-hnear evolution is the Press & Schechter (PS, 
Press & Schechter 1974) approach. It is based on the spherical collapse model, 
that Gunn & Gott (1972), Gott & Rees (1975), Peebles (1980) debated within 
the frame of pure CDM models, and Lahav et al (1991), Eke et al (1996), 
Brian & Norman (1998) and others generalized to the case of ACDM. In 
spite of its approximation, such model, inserted in PS formulation, has been 
found to approach simulation results (see, e.g., Lacey & Cole 1993, 1994). 
Recent improvements of the method (Sheth & Tormen 1999, 2002, Sheth, Mo 
& Tormen, 2001, see also Jenkins et al 2001), allowing a better approximation, 
involve some more parameters and their use seems unnecessary when aiming 
just to compare different cosmological models. 



2 From linear theory to non linecir predictions 

In order to apply the PS technique we first need to determine the ampli- 
tude 5c, in the linear theory, of a spherical fluctuation that would achieve full 
recoUapse at a given redshift Zcoi- Real fluctuations, after achieving maximum 
expansion, turn around and begin to recontract. However, contraction requires 
that potential energy, turned into internal kinetic energy, is radiated soon. If 
radiation is negligible, contraction is delayed and, at Zcoi, virial equilibrium 
is attained. In a standard CDM model, Sc is constant and holds ~ 1.68 (see, 
e.g.. Coles & Lucchin 1995). In RP or SUGRA models, Sc depends upon the 
cosmological parameters and the redshift Zcoi- In Fig.s 1 and 2 we report such 
dependence. As a by-product of such computation, we obtain the density con- 
trast (Ac) of a spherical fluctuation when it is fully virialized. In Fig. 3 we 
report the dependence of Ac on ai z — 0, for a number of DE models. 
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Fig. 1. The dependence of 5c on the matter density parameter is shown, for 4 
RP (A/GeV = 10^, 10^, 10^ and 10^) and 2 SUGRA models (A/GeV = 10^ and 
10^). A values increase from top to bottom curves. 

Then, in Fig. 4, we report how Ac depends on z for some DE models, keeping 
Qm = 0.3 and h = 0.7. 

The value of Sc is to be used in the expression 

M 

f{u)udlogu = —nc{M)MdlogM ; (2) 

Pm 

(here is background matter density, the bias factor u = Sc/(Jm, ctm being 
the rms density fluctuation on the length scale corresponding to the mass M) , 
yielding the differential mass function ndM) once the distribution on bias 
is given. Here, as usual, we assume a Gaussian /(z/). Eq. (2) shall then be 
integrated, to obtain nc{> M). Let us outline that, in order to compute ctm, 
we need to know the transfered spectrum P{k)] let us also recall the relation 
k = 27r(47rp/3M)^/^ between wavenumber and mass. 

The first finding of this analysis is that the difference between the mass func- 
tions for RP or SUGRA models and the mass function of a AGDM model, at 
2; = 0, is smaller than the difference due to a shift by 0.05 of the primeval 
spectral index. The intrinsic noise of data is then a serious obstacle to any 
attempt to determine the nature of DE using cumulative mass functions at 
^ = 0. 

On the contrary, we found a clear imprint of the nature of DE in the evolution 
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Fig. 2. The dependence of dc on the redshift z is shown for $7^ = 0.3 and h = 0.7. 
A values and curve setting are the same as in Fig. 1. 
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Fig. 3. The density contrast Ac of a virialized sphere is shown for 4 RP and 2 
SUGRA models as a function of 0,rm for h = 0.7, at z = 0. A values are the same 
as in Fig. 1, but here they decrease from top to bottom curves. 

of the number N of clusters (above a suitable mass M, in a box of side s). 
In Fig. 5, we give (for M > 6.9 • lO^^/i-^M© and s = 100 h'^Mpc) as a 
function of redshift, normalized to an identical number of clusters at 2; = 0, 
for all models {N(z = 0) = 0.13). In Fig. 6 we give the ratio between the 
number of clusters expected in various models and the number expected in 
an open CDM model with the same value of Qm- The mass here is selected 
so to correspond to a cluster with Abell radius in a standard CDM model. 
A similar plot, for a slightly smaller mass, was given by Bahcall, Fan & Cen 
(1997), for standard CDM, ACDM and OCDM only (for a recent analysis of 
the constraints that cluster number counts set to the cosmological model, see, 
e.g.. Holder, Haiman & Mohr, 2001). 



3 Discussion 

The above mass functions and evolution are calculated using a transfer func- 
tion obtained from a generalization of the public program CMBFAST to cos- 
mologies with DE given by the potentials (1). Initial conditions were set ac- 
cording to the tracker solution in radiation dominated era (Steinhardt, Wang 
& Zlatev 1999; Zlatev, Wang, & Steinhardt 1999, Brax, Martin & Riazuelo 
2000). DE fluctuations, taken into account in primeval fluctuation evolution, 
are damped soon after a scale enters the horizon and, therefore, are not im- 
portant in further fluctuation evolution. Accordingly, in our problem, their 
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Fig. 4. The density contrast of a virialized sphere is shown for 4 RP and 2 
SUGRA models as a function of z. For the sake of comparison, the constant value 
of Ac in standard CDM is also reported. A values and curve setting are the same 
as in Fig. 3. 
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relevance amounts to causing some modification of the transfer function. (On 
the contrary, they cause major modifications of CMB anisotropy and polar- 
ization.) 

Fig. 7 describes the pattern followed to evaluate 5c and for Zcoi — 0; for 
Zcoi > the procedure is similar. It must be however outlined that Fig. 2 does 
not describe a time evolution, but is worked out considering different time 
evolutions of the kind described in Fig. 7. We obtain the linear growth factor 
by using a simplified program, suitable to treat fiuctuation evolution after 
recombination, which includes the Friedman equation, the equation of motion 
of the scalar field and the Jeans' equation for the linear density fiuctuation 5: 

U2lS-I^S^0 (3) 

(dots yield derivatives in respect to time. Ho is the present value of the Hubble 
constant). Such system provides the lower curves in Fig. 7; for a fiat model 
without DE (SCDM), the curve is a straight line; together with it, in the 
figure, we show the behavior for a RP model. Other models, with different 
DE contents, have a similar behavior. Such linear results can be multiplied by 
an arbitrary factor and, in the figure, they are normalized so that the initial 
fiuctuation amplitude allows full recoUapse, at Zcoi, of a spherically symmetric 
fiuctuation. This is also to be done in actual calculations to evaluate 5c- The 
evolution of the contrasts A = Pm/Pm, between inner and average matter 
densities, is then shown by the upper curves. However, instead of reporting 
their diverging at Zcoi, we interrupt their growth when virialization density 
contrast is attained; the residual slight growth of A is due to cosmic expansion. 
Notice that the density contrasts Ac, shown in Fig.s 3 and 4, are il^A. 

More in detail, within such spherical lump, at an initial redshift Zi, the energy 
density exceeds average by a factor 1 + 5i — Ai, tuned so cause full recoUapse 
at an assigned redshift Zcoi — i/cicoi — 1- The radius R of the lump starts 
from Ri] its evolution is then computed by taking into account, besides of the 
varying energy density (pm) of the baryon and dark matter within R, also 
the energy density and pressure of DE {pde and Pde, respectively), within R 
itself. In models involving just a cosmological constant, it is —poE = Pde ~ 
const., during the whole process (see, e.g., Lahav et al 1991). In the present 
case, instead, pde, Pde and their ratio vary in time, simply because of their 
overall evolution with a. At turn-around the radius will be Rf^. Afterwards, 
R decreases and should vanish at Zcoi- The ratio between the energy density 
inside R and average yields A. 

For the sake of example, in Fig. 8 we show R(t) for two RP potentials (A/GeV = 
10^ and 10^) and ACDM. SUGRA potentials yield curves intermediate between 
RP and ACDM models and just slightly dependent on A. A careful inspection 
of the figures shows that, at variance form SCDM and ACDM, RP models 
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Fig. 5. The number density of clusters is shown for standard CDM (lowest curve), 
OCDM (long-short dashed upper curve), ACDM (solid intermediate curve), 2 RP 
(long dashed curves) and 2 SUGRA (dotted curves) models. The values of A/GeV 
are 10^, 10^ for RP and 10®, 10^*^ for SUGRA, starting from the lower curves. 




Fig. 6. The number densities of the previous plot normalized to the values in OCDM 
are shown. 
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Fig. 7. Linear and non-linear evolution of density fluctuations. A models with DE 
due to a self-interacting scalar field is compared with SCDM. Similar plots can be 
given when full recollapse is expected when a < 1. The density contrasts Ac, shown 
in Fig.s 3 and 4, are ^m^- 

have a slight asymmetry between expansion and recontr action, due to the 
evolution of DE density. In Fig. 8, the virialization radius Rvir, from which Ac 
is computed, is also indicated. In the SCDM case, energy conservation yields 
Ryir = Rtaf^- In the presence of DE, requiring virial equilibrium and energy 
conservation leads to the cubic equation 



2 + yjata) 



X + } , = . (2) 
2r;(ac„^) ^ ^ 



Here x — Rvir I Rta and 



^(^) = 2^.(a)(l + A.) [-R-) UJ 

is the matter density parameter when the scale factor is a; accordingly, 
^2^(1) = flm)- Solving such equation yields x values slightly below 0.5 . 



4 Conclusions 



Observable effects of the nature of DE have been considered by various au- 
thors. For instance, Cooray & Huterer (1999) discussed the relation between 
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Fig. 8. Time dependence of the radius of a density fluctuation in various models; 
horizontal lines, at the right, indicate the Rmr/Ri attained at to- Similar plots can 
be given when full recollapse is expected at t <to. 
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Fig. 9. Rcdshift dependence of w for 4 RP (A/GeV = 10^, 10^, 10^ and 10^) and 
2 SUGRA models (A/GeV = 10^ and 10*^). A values decrease from top to bottom 
curves. 

DE nature and gravitational lensing. Previous work on the value of 6c was 
made by Steinhardt, Wang & Zlatev (1999), although explicit outputs were 
not given. More recently, Lokas (2002) considered the behavior of Ac and the 
mass function in DE models with constant w — —p/p- This approximated 
treatment has been pursued in a number of recent papers and eases compu- 
tations. Let us however remind that the value of to, when dynamical DE is 
considered, varies significantly in the relevant period. In Fig. 9 we report its 
variations between z — Q and 10, when structures form. Notice that the rate 
of w variation is highest in the most critical redshift interval, between z — Q 
and 1-2. In RP models such variation is ~ 20 %. In SUGRA models it is even 
greater, as w passes from values ~ 0.8, at 2; = 0, up to values ~ 0.3-0.4 
in the above narrow z interval. Taking constant therefore, is a dangerous 
approximation, whose reliability ought to be carefully inspected, in different 
problems. 

The main results of this work are that: (i) the shape of the mass function of 
clusters, at 2; = 0, is only mildly dependent on DE nature. On the contrary, (ii) 
the cluster evolution depends on the nature of DE in a significant way. More 
in detail, models with RP potentials closely approach the evolution expected 
in open CDM models. Only for A values as low as ~ 10^ GcV, the expected 
behavior in a RP model is appreciably different from an open CDM with the 
same ilm- On the contrary, the evolution of SUGRA models is intermediate 
between open CDM and ACDM models. 
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Cluster data available within a few years were thought to be able to discrimi- 
nate between open CDM and ACDM, on the basis of the redshift dependence 
of cluster abundance. If independent data confirm that we live in a spatially 
fiat world, finding an evolution closer to open CDM than to ACDM will pro- 
vide a precise information on the nature of DE. 

As a by-product of the analysis leading to these conclusions, we also found 
the dependence of the virialization density contrast Ac on flm, on the nature 
of DE, and on redshift z. Such Ac is to be used in various applications, e.g. 
to build SO algorithms able to find clusters in n-body simulations of models 
with DE. 
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